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ABSTRACT 


The  statistical  theory  of  radio  direction  finding  as  it  applies  to  three-dimensional  geometries  is 
presented  in  this  report.  This  theory  is  an  extension  of  the  two-dimensional  theory  presented  by  R.  G. 
Stansfield  in  1947.  The  theory  has  application  in  situations  where  airborne  targets  are  maneuvering  in 
three  dimensions  at  high  speed;  and,  especially  when  the  platforms  with  the  direction  finding  equip¬ 
ment  are  also  aircraft. 


:  ?  IV3'™"'  '&•?'  '^mt<^!f!fi^pvr~ 


IV',  -  s4  W^>4  •'.'■•V,  '>  *  —  «-. 


THE  JOHNS  HOPKINS  UN^ERSITY 

APPUEO  PHYSICS  LABORATORY 
laurel  Maryland 


Section 


CONTENTS 


List  of  Illustrations . . .,,,„,, . . . .,, ,:. ., . 


Introduction., 


Stansfield’s  Theory 


2.1  The  Two-Dimensional  DF  Encounter .,,,, . .,,, 

2.2  The  Two-Dimensional  DF  Analytical  Results .,,, .. . . .,,. 


The  Three-Dimensional  Extension  of  Stansfield’s  Theory 

3.1  The  Three-Dimensional  DF  Encounter.;., . . . 


3.1.1  Initial  Concepts  . . . 

3.1.2  The  Construction  of  p,  and  </,.... . 

3.1.3  The  Fundamental  Statistical  Characterization  of  the 

Three-Dimensional  DF  Theory 

References..:..,...,, . . . ., .,, . . ., .,, , 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 
laurel.  Maryland 


ILLUSTRATIONS 


Figure  Page 

2-1  Two-Dimensional  Radio  DF  Geometry  and  Para  me  ters . ... ... ............ ... ., ..., ., . .... . ......  2 

2-2  A  Magnified  View  of  the  Coordinate  Origin  to  Determine  q,  ....... . 4 

2-3  The  Coordinate  System  ,x«jo  Rotated  About  O  Relative  to  the  System  . .  8 

2- 4  A  Contour  of  Equal  Likelihood  Bounding  the  Region  .4^.,,,,,. . . 10 

3- 1  Three-Dimensional  Radio  DF  Geometry  and  Parameters..:.:.:., ,,,,, 14 

3-2  A  Magnified  View  of  the  Coordinate  Origin  to  Determine  q, 

When  cjj  and  p,  are  Orthogonal  to  /3; . 14 

3-3  A  Magnified  View  of  the  Coordinate  Origin  to  Determine^  q, 

When  qj  is  Orthogonal  to  Dj+7 and  p,  is  Orthogonal  to  D, 15 

3-4  The  («,  fi,y) ,  Coordinate  Axes  Relative  to  the  (.v,v,r)  Coordinate  Axes.. .... ... . . .... . 20 

3-5  The  Resolution  p,  in  the  (fS-y)j  Plane  of  the  (a.fi.y) ,  Coordinate  System ...... 20 

3-6  Graphical  Inversion  of  P{x'q  r'o.zV*  . ....... . . ., .,,, 36 


-  vii  - 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL  MARYLAND 


SECTION  1 
INTRODUCTION 


It  is  the  primary  objective  of  this  report  to  present  the  principal  features  and  results  of  a  statistical 
theory  descriptive  of  the  process  of  radio  direction  finding.  The  theory  presented  here  is  based  upon 
that  developed  by  R.G.  Stansfield  as  documented  in  his  paper  of  1947  (Reference  1).  The  unique 
feature  of  the  theory  presented  herein  is  its  applicability  to  three-dimensional  geometries,  as  opposed  to 
the  two-dimensional  nature  of  the  Stansfield  theory;  in  this  sense,  then,  the  radio  direction  finding 
theory  presented  within  this  report  is  an  extension  of  Stansfield's  theory. 

Interest  in  pursuing  an  extension  of  Stansfield's  theory  to  three  dimensions  was  aroused  as  the 
result  of  a  discussion  with  a  colleague  regarding  the  problem  of  passively  determining  the  position  of  a 
moving,  radiating  target.  It  became  obvious  that  >t  would,  in  general,  be  necessary  to  employ  several 
dispersed  direction  finding  sites  capable  of  simultaneous  bearing-line  measurements. 

Conceptually,  this  circumstance  leads  to  no  difficulties  until  one  begins  to  consider  applying 
Stansfield's  theory  to  provide  analytical  characterization  of  such  a  direction  finding  scheme.  Indeed,  if 
one  contemplates  the  situation  wherein  the  target  is  a  maneuvering,  high-speed  aircraft,  and  the  direc¬ 
tion  finding  sites  consist  of  several  "picket"  or  early  warning  aircraft  with  direction  finding  equipment 
onboard,  it  would  be  very  difficult  if  not  impossible  to  control  the  picket  aircraft  so  that  they  and  the 
target  aircraft  stayed  in  the  same  geometric  plane  during  the  period  of  encounter.  This  would  never  be 
done  operationally;  but  one  would  have  to  assume  this  geometry  to  apply  Stansfield's  theory  directly  to 
any  analytical  evaluation  of  the  encounter.  This  is  not  meant  to  suggest  that  Stansfield's  theory  has  no 
useful  applications;  it  has,  in  fact,  been  the  cornerstone  of  direction  finding  analysis  for  the  last  33 
years.  In  each  instance  where  it  has  been  applicable,  however,  the  physical  circumstance  under  investi¬ 
gation  has  allowed  valid  application  of  the  two-dimensionai  theory  (e.g.,  see  References  2  and  3).  It  is 
still  true  today  that  the  bulk  of  the  hardware  devoted  exclusively  to  radio  direction  finding  applications 
can  measure  bearing  lines  in  a  single  plane  only.  The  point  to  be  made  here  is  that  there  are  a  number 
of  circumstances  in  which  the  ability  to  perform  three-dimensional  direction  finding  measurements 
would  be  extiemely  useful,  and,  to  that  end,  a  three-dimensional  extension  of  Stansfield's  two- 
dimensional  direction  finding  theory  is  necessary  in  order  to  provide  a  suitable  analytical  tool  for 
evaluating  the  effectiveness  of  current  or  prospective  three-dimensional  direction  finding  systems. 

Circumspect  review  of  this  document  by  persons  cognizant  of  or  interested  in  direction  finding 
applications  and  procedures  is  deemed  essential  and  is  heartily  invited.  Comments,  suggestions,  and 
criticism  are  encouraged  not  only  to  support  the  evolution  of  a  more  intelligible,  useful  report,  but  also 
to  help  m  assessing  the  basic  wer’h  of  and/or  need  for  what  is  presented  here.  Indeed,  there  is  some 
doubt  as  to  the  uniqueness  of  this  material,  i.e.,  even  though  a  limited  rev  iew  of  some  current  litera¬ 
ture  (e.g.,  see  Reference  4)  has  failed  to  uncover  or  suggest  any  other  three-dimensional  extension  of 
Stansfield's  direction  finding  theory,  one  cannot  help  but  feel  that  this  should  have  (and  certainly  could 
have)  been  accomplished  prior  to  this  date. 

The  remainder  of  this  memorandum  will  present  first  a  brief  discussion  of  Stansfield's  two- 
dimensional  radio  direction  finding  theory,  and  then  present  a  discussion  of  the  three-dimensional  radio 
direction  finding  theory.  A  precedent  discussion  of  Stansfield's  theory  will  obviate,  at  least  initially,  the 
acquisition  oi  Stansfield's  original  paper  by  the  reader  and  will  provide  a  convenient  point  of  departure 
from  which  to  initiate  a  discussion  of  the  three-dimensional  theoiy.  Of  even  more  importance,  how¬ 
ever,  is  the  sense  of  correspondence  and/or  distinction  to  be  gained  by  the  reader  when  comparing 
Stansfield's  two-dimensional  results  with  those  of  the  thiee-dimensional  theory;  e.g..  as  a  minimum 
requirement  for  credibility,  cue  would  expect  the  results  of  the  thiee-dimensional  theory  to  coincide 
with  those  of  Stansfield’s  theory  for  planar  geometries  with  null  elevation  error,  and  to  appreciate  such 
coincidence,  the  reader  must  be  acquainted  with  Stansfield's  theory. 

In  what  follows,  the  terms  "direction  finding,”  or  "direction  finder"  will  be  abbreviated  as  T)F," 
where  the  meaning  should  be  clear  from  the  context  in  which  the  abbreviation  is  used.  Likewise  "BL" 
and  "BA"  will  be  used  to  represent  the  words  "bearing  line"  and  "bearing  angle,"  respectively. 
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SECTION  2 

STANSFIELD’S  THEORY 
2.1  THE  TWO-DIMENSIONAL  DF  ENCOUNTER 

Figure  2-1  illustrates  the  general  distribution  of  the  participants  and  identifies  the  parameters  used 
to  characterize  the  two-dimensional  DF  encounter.  The  target  transmitter  whose  location  is  to  be 
determined  and  all  of  the  DF  sites  attempting  to  measure  the  angles  of  BLs  to  the  transmitter  are  Con¬ 
sidered  to  be  in  the  same  geometrical  plane.  Points  in  this  "encounter  plane"  are  labeled  via  a  fixed 
cartesian  coordinate  system  with  an  origin  Oat  the  position  of  the  target  transmitter. 

As  mentioned,  each  DF  site  employs  some  form  of  directional  antenna  system  in  order  to  derive 
the  angle  of  a  BL  from  the  DF  site  toward  the  transmitter  location.  The  situation  forjhe  j’th  DF  site  is 
shown  in  detail  in  Figure  2-1.  The  BL  to  the  actual  transmitter  location  is  the  line  JO  with  a  BA  of  0, 
and  length  D,._However,  due  to  instrumental,  propagationai,  and  operator  errors,  the  measured  BL  is 
along  the  line  ./Pwith  a  BA  of  +  t//,.  where  ip,  is  the  BA  error. 


Figure  2-1  Two-Dimensional  Radio  DF  Geometry  and  Parameters 


A  fundamental  assumption  made  by  Stansfield  is  that  the  BA  error  t b,  is  the  value  of  a  normally  distri¬ 
buted  random  variable  'J'^,  where 


P(<!>,  a  V.  ^  i//;  +  </./»,)  =  pi# ,) ciifr ,  =  — exp 

<r,v 2n)  ■  [Ivj 


dilij 


(1) 
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A  similar  assumption  will  be  made  for  the  azimuth  and  elevation  BA  errors  when  three-dimensional 
geometries  are  considered. 

Rather  than  use  i|/(  as  the  basic  indication  of  the  amount  by  which  the  j'th  DF  kite's  BL  measure¬ 
ment  is  in  error,  the  length  p ,  of  the  line  PO  is  used.  The_  line  PO  is  constructed  from  the 
transmitter  location  O  perpendicular  to  the  measured  BL  along  JP.  The  length  p ,  is  called  Ihe  bearing 
error.  Another  fundamental  assumption  made  by  Stansfield  is  that  the  standard  deviation  of  the  BA 
errors  or,  is  so  small  that  the  region  of  uncertainty  determined  by  the  entire  set  of  measured  BLs  is 
small  in  comparison  with  the  distances  Dr  In  this  case,  p ,  is  related  to  >h,  by  the  approximate  relation 


P,  =  /->,*/», 

This  means  that  p ,  is  the  value  of  a  random  variable  Pr  where 


(2) 


1  \ 

p 

p.  dP, 

\p,  «:  p,  +  dp,\  =  P 

D, 

^  D,+  D, 

=  /»</>,)<//>,  =  — rr;  exp 
<r  Alir)  ■ 


-rr 

2a/- 


dp, 


(3) 


and  where 


<rj  =  <r,D, 


(4) 


One  of  the  n  DF  sites  of  Figure  2-1  is  also  a  DF  operations  center  to  which  all  of  the  other  DF 
sites  send  reports  giving  their  measured  BAs  to  the  target  transmitter.  Having  received  the  complete 
set  of  BA  data,  the  personnel  of  the  DF  operations  center  must  somehow  use  that  data  to  generate  an 
estimate  of  the  location  (the  coordinates)  of  the  transmitter.  Suppose  that  the  operations  center  per¬ 
sonnel  hypothesize  the  transmitter’s  location  to  be  at  the  point  S.  Relative  to  the  fixed  reference  sys¬ 
tem  with  origin  at  O,  the  line  OS  defines  a  position  vector  7  =  xi\  +  vh,,  so  that  the  coordinates  of  the 
point  5  relative  to  the  reference  system  are  (.v,.i).  It  must  be  emphasized  here  that  the  DF  operations 
center  personnel  have  no  knowledge  of  the  coordinates  (.v..i )  ol  the  point  5  relative  to  the  reference 
system  at  0,  for  if  they  did,  they  could  maneuver  their  hypothetical  point  5  so  as  to  diminish  .v  and  ,t 
and  locate  the  transmnter  with  an  arbitrarily  small  error.  The  coordinate  reference  system  at  0  exists 
for  analytical  convenience  and  is  not  intended  to  represent  any  operational  system  of  measurement. 

However,  when  the  BL  for  any  DF  sight  is  plotted,  the  DF  operations  center  personnel  can  meas¬ 
ure  the_  'perceived  error."  In  Figure  2-1,  the  measured  BA  of  the  j’th  DF  site  yields  a  BL  that  lies 
along  JQ  and  the  perceived  error  is  the  length  of  the  perpendicular  line  SQ,  written  as  qr  The  per¬ 
ceived  error  qr  in  contradistinction  with  the  actual  error  p„  is  the  bearing  error  attributed  to  the  j’th 
DF  site  by  the  personnel  of  the  DF  operations  center  as  a  result  of  their  hypothesis  that  the  transmitter 
is  at  point  S.  Of  course,  q ,  and  p,  are  further  distinguished  by  the  fact  that  q,  is  a  quantity  known  to 
the  personnel  of  the  DF  operations  center,  while  p,  is  not  known  to  them.  Indeed,  the  DF  operations 
center  personnel  can  employ  knowledge  of  the  set  of  perceived  errors  {</;},  A  =  1,  2,  3  n  (one  for 
each  DF  site)  to  select  a  "best  estimate"  or  "fix"  for  the  hypothesized  location  S of  the  transmitter. 


The  set  of  perceived  errors  UJ  accumulated  at  the  DF  operations  center  is  considcied  to  be  a 
sample  of  n  independent,  normally  distributed  random  variables  Q\.Q:  Q„.  In  this  case,  the  proba¬ 

bility  of  observing  any  specific  set  of  perceived  errors  \qj  is  given  by: 
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The  best  estimate  for  the  location  of  the  transmitter  is  now  taken  10  be  the  point  S whose  associated  set 
of  perceived  errors  has  the  highest  likelihood  of  being  observed.  As  mentioned  by  Stansfield,  this 
is  an  application  of  the  statistical  sampling  principle  known  as  the  theorem  of  maximum  likelihood. 
The  expression  of  equation  (5)  will  be  maximized  by  that  set  of  perceived  errors  [qk\  that  minimizes 
the  argument  A  of  the  exponential,  where  ■ 


A  - 


(6) 


If  we  now  define  ea  = 


(the  vertical  bars  signify  the  absolute  value  of  the  ratio  between  them) 


as  the  relative  bearing  error  of  the  k’th  DF  site,  the  quantity  A  is  seen  to  be  the  sum  of  the  squares  of 
the  relative  bearing  errors  over  all  the  DF  sites.  Stansfield’s  maximum  likelihood  criterion  for  choosing 
the  fix  point  .Vof  the  transmitter  is  seen  to  be  equivalent  to  selecting  the  point  S' whose  associated  per¬ 
ceived  errors  UJ  satisfy  a  ’  least  squares"  criterion  as  applied  to  the  relative  errors  ek. 


In  Figure  2-2,  a  unit  vector  q,  (horizontal  bars  above  single  letters  signify  a  vector  quantity,  and  a 
carat  above  a  letter  indicates  that  the  quantity  being  described  is  a  unit  vector)  is  defined  along  the 
direction  of  the  line  QS  and  is  given  approximately  by 


cij  =  sin  -  cos  0,u,  (7) 

and  reference  to  the  figure  yields  for  q,  (the  length  of  QS )  the  expression 

<7,  =  P,  +  ~r  ■  q,  (8) 

The  "dot"  between  r  and  q ,  indicates  a  scalar  product  of  these  vectors.  Employing  equations  (7)  and 
(8),  equation  (6)  can  be  rewritten  as 


PO-75 


Figure  2-2  A  Magnified  View  of  the  Coordinate  Origin  to  Determine 
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n  ~ 

E(p,  +  a  sin  0,  -  y  cos  tt,r 

Ur,Dy 


(9) 


The  sum  of  the  squares  of  the  relative  errors  is  now  expressed  in  terms  ot  the  actual  bearing  errors  n, 
and  the  coordinates  (a.i  )  of  the  hypothesized  transmitter  location  S.  This  is  a  critical  result  that 
allowed  StansfieJ  to  find  closed  form  expressions  for  the  fix  coordinates  l.x.i  )  for  any  single  attempt  at 
locating  the  transmitter  when  the  actual  bearing  errors  are  given  and  to  describe  the  statistical  features 
of  the  pair  (a.i)  when  the  actual  bearing  errors  are  characterized  as  random  variables. 


2.2  THE  TWO-DIMENSIONAL  DF  ANALYTICAL  RESULTS 


If  one  uses  equation  (9)  to  evaluate 


0.4 


O.v 


and 


04 


by 


,  then  the  equations 


bA 

04 

bx 

=  0  and 

9.v 

=  0 


constitute  two  equations  in  the  two  unknowns  .vand  r  which  can  be  solved  for  a  and  .i  Solving  these 
equations  one  finds 


where 


1 

(\/Lt  —  I  >2) 


i ’  cos  ft,  -  ix  sin  ft. 


(10) 


1 

(\(JL  —  l'?) 


P, 


( <r  ,D,)~ 


~(\  cos  ft,  -  p  sin  ft,) 


(11) 


\  = 


sin*  ft , 
(ir  ,D,)2 


_  ^  cos2  ft, 
(f r  ,D,)- 


^  sin  ft,  cos  ft, 

,,  —  ^  — — - — — 

Ur  ,D,)2 


(12) 

(13) 

(14) 


Although  not  represented  explicitly,  the  summations  are  over  the  index  ./for  j=  1,  2  ....  n;  i.e,  the  sum¬ 
mations  contain  one  term  associated  with  each  DF  site 

Equations  (10)  and  (11)  are  Stansfield's  equations  for  the  fix  coordinates  that  would  be  deter¬ 
mined  by  the  personnel  of  the  DF  operations  center  for  any  single  location  attempt;  wherein  the  actual 
bearing  errors  are  the  set  {/>,!,/  =1,2  n.  Other  than  the  set  of  actual  bearing  errors,  the  values  for 
the  fix  coordinates  depend  upon  the  geometrical  distribution  of  the  DF  sites  relative  to  the  actual 
transmitter  location  (represented  by  the  parameter  sets  {ft,}  and  {D,})  and  the  instrumental,  propaga- 
tional,  and  operator  error  characteristics  (represented  by  the  parameter  set  {cr,)). 

Notice  that  if  we  do  not  wish  to  or  cannot  describe  the  distribution  of  the  bearing  errors  as  nor¬ 
mally  distributed  about  the  actual  target  position,  we  can  use  a  more  general  form  of  equation  (6)  given 
by 


A  =  y  q 

*  \VA  / 


(15) 


where  '  Q k  2  is  the  mean  squared  value  of  the  random  variable  ft  which  characterizes  the  bearing  error 
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of  the  k’th  DF  site.  If  Qk  is  other  than  normally  distributed  about  the  actual  target  position,  then,  in 
general,  (Qk2)  will  not  be  equal  to  J  .  In  this  case,  equations  (10)  through  (14)  are  still  useful 
when  (  )2  is  replaced  by  \Qk')- 

If  one  rewrites  equations  (10)  and  (11)  as 


1 

y  Pj  cos  Vj 

-  M 

y  P,  sin  0, 

(A/i  -  v2) 

V 

( <r,D ,)2 

( <TjD ,)2 

1 

V 

y  p,  cos 

V  P,  sin  0, 

1 

4 

-< 

A 

(rr;D;)2 

—  V 

(< TjD ,)2 

(16) 


(17) 


one  sees  that  these  equations  take  the  form 


where 


.  _  _ v _ _ /i 

( A/x  -  v2)  (A/i  —  v2) 


A  i' 

y  =  - yr  u  -  - - rr  t 

(A/i  -  i'2)  (A/i  —  jj2) 


y  p,  cos  u, 

"  ^  ( cr,D,)2~ 


y  P,  sin  0, 
'  ^  (<r,D,)2 


(18) 

(19) 


(20) 


(21) 


Since  pt  is  the  value  of  a  random  variable  Pr  it  is  clear  that  ii  and  t  are  the  values  of  random  variables 
we  will  designate  as  Y  and  <9,  respectively,  and  thus  x  and  y  are  the  values  of  random  values  to  be 
designated  as  X  and  K  respectively.  Equations  (18)  and  (19)  mean  that  our  random  variables  satisfy 
similar  equations;  i.e.. 


X  = 


Y  - 


(A  /i  —  i'2)  (A/i  —  v2) 


0 


(22) 


)  = 


Y  - 


(A/i  —  i'2)  (A/i  —  v2) 


0 


(23) 


Recall  that,  consistent  with  Stansfield's  assumptions,  the  random  variables  {  P,  }  are  all  distributed  in  a 
similar  fashion  with  zero  means.  Thus,  the  central  limit  theorem  can  be  invoked  to  contend  that  the 
random  variables  Y  and  0  are  normally  distributed  with  zero  means  (Reference  5).  This  being  the 
case,  we  also  know  that  X  and  Y  are  normally  distributed  (reference  6). 

If  Y  and  0  are  normally  distributed,  then  their  joint  density  will  be  given,  in  general,  by  (Refer¬ 
ences  5,  6,  and  7) 
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- 1  ir  2  Rut  ,  i 2 

exp - —  - r - +  - 7 

|  2(1  -  R2)  it 1 2  a  2 

PiU,,)  ~  2ir  rr „  a,  (1  —  RrY/: 


where  <r„  and  <r,  are  the  standard  deviations  of  Y  and  0,  respectively,  and  R  is  the  correlation 
coefficient  of  Y  and  (-).  Evaluation  of  <r„.  a-,,  and  R  give 


'/• 

<r„  —  n  ■ 


(T,  —  A’ 


If  equations  (18)  and  (19)  are  inverted,  they  yield 


u  =  fxy  —  vx 


t  =  w  —  kx 


and  substitution  of  equations  (25)  through  (29)  into  (24)  gives 


p'(x.y)  = 


exp  -y(A.r  -  2uxy  +  /xy)| 
2n(kn  —  u2)V: 


When  changing  coordinates  from  the  (u,i)  space  to  the  (x,y)  space,  we  have  in  general  that 

p(x,v)  =  p'(x,y)J(u.l/x,y ) 

where  J(u,t/x,y)  is  the  Jacobian  determinant  for  the  transformation..  In  this  case, 

J(u.i/x,y)  =  Aft  —  v2 

so  with  (31)  and  (32),  equation  (30)  becomes 


p(x,y)  =  -  -  exp  -—(A.v2  -  2i>xv  +  fi\2) 

2ir  2 


Following  equations  (iO)  and  (11),  equation  (33)  constitutes  the  second  significant  analytical 
result  of  Stansfield’s  two-dimensional  radio  DF  theory.  Equation  (33)  is  the  joint  probability  density  of 
the  random  variables  A' and  Y,  the  coordinates  of  the  DF  position  fix  for  the  target  transmitter.  When 
multiplied  by  the  differential  area  of  a  neighborhood  about  the  point  (.x;.y),  it  determines  the  likelihood 
that  the  DF  fix  resulting  from  any  given  attempt  to  locate  the  target  transmitter  will  lie  within  that 
differential  neighborhood  of  the  point  (x,y ).  The  quadratic  nature  of  the  argument  of  the  exponential 
and  the  fact  that  A  and  ft  are  greater  than  or  equal  to  zero  indicates  that  the  contours  of  constant  likeli¬ 
hood  are  ellipses  in  the  encounter  plane  centered  on  the  actual  target  position  (the  point  0  in  Figures 
2-1  and  2-2).  It  is,  of  course,  possible  to  specify  a  system  of  coordinates  rotated  about  Dby  the  angle 
<f>  relative  to  the  system  x-y  in  terms  of  which  these  elliptical  contours  can  be  expressed  as  a  simple 
sum  of  squares  of  the  coordinates  (A:0,y0).  If  in  Figure  2-3  we  suppose  that  the  locus  of  points  which 
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constitute  a  given  elliptical  contour  of  constant  likelihood  can  be  expressed  in  canonical  form  as 


*L.  +  K.  =  k2 

a2  b2 


where  k2  is  a  positive  constant,  then  use  of  the  coordinate  transformation  between  the  a-o^o  and  x  y 
coordinate  systems  with  equation  (34)  yields  an  expression  for  the  locus  as  expressed  in  the  x  y  coordi¬ 
nate  system 


Figure  2-3  The  Coordinate  System  x0y0  Rotated  About  O  Relative  to  the  System  xy 


cos*  <b  ,  sin2  cb  i  ,  .  ,  .1  1 

- —  +  — —  r  +  2  sin  <b  cos  <b  —r-  -pr  XV 

a  bl  a 2  bl 


sin2  <j>  ,  cos' 


+  v2  =  c2 


Comparison  with  the  argument  of  the  exponential  in  equation  (33)  yields 

_  cos2  <b  ,  sin2  & 


v  =  —  sin  <j)  cos  <t>  I  -4 - M 

l  a  b  ) 

sin2  ,  cos2  <t> 

/x  =  — r —  +  — 

a i  b 


From  equations  (36)  through  (38)  we  find 


-9- 


r*tt 


I 

I 

1 
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A  family  of  concentric  ellipses  about  the  point  O  to  which  the  ellipse  of  equation  (50)  belongs  is 
delined  by 

^4-+^-=  K2.  0  <  K  <oo  (51) 

a2  b2 

The  probability  that  the  fix  point  (  .y0..Vo  )  will  lie  between  the  elliptical  contours  defined  by  K  and  K  + 
dK  is,  by  using  equation  (47), 


J 

ATa + (IK 


p(x0,y0)  dx0dyQ  =  K 


exp 


K± 

2 


JK 


Thus,  the  desired  probability  can  be  evaluated  by 


K2 

[■Vo  ,.v 

0  =  j  K 

'  0 

exp 

2 

(52). 


(53) 


so 


r  , 

,  ,  j 

k2 

[.vo  > 

o  —  1  - 

exp 

2 

(54) 


Since  the  coordinates  of  the  locus  of  points  defining  the  boundary  of  the  region  Ak  satisfy  equation 
(50),  equation  (54)  can  be  rewritten  as 


[  ■  A  , 

1 

-vo2 

7 

,  To' 

(■VoJ’o  €  j  =  1  - 

exp 

2 

■> 

a  ~ 

(55) 


Subtracting  1  from  both  sides  and  taking  the  natural  logarithm  (In)  of  both  sides  gives  another  of 
Stansfhld's  results 


I 

I 

1 


1 

£T 


+ 


(56) 


Equation  (56)  expresses  a  convenient  relationship  between  the  locus  of  points  bounding  an  elliptical 
region  Ak  and  the  probability  that  the  coordinates  of  a  given  DF  fix  will  fall  inside  Ak.  It  is  often  desir¬ 
able  to  specify  a  required  value  of  P(x0',va'eAk)  and  to  use  (56)  to  determine  the  resultant  dimensions 
of  the  region  Ak.  Indeed,  the  lengths  of  the  semimajor  and  semiminor  axes  (or  vice  versa  if  b  >  a) 
are 


UoJ 


1 

l-P(x0.y0€Ak) 


ITO/h  1 


1 

l-P{x0,y0eAk) 


-  II  - 


(57) 


(58) 
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and  the  a  ea  of  the  region  An  is 


=  nab  In 


\-P(x'0,yoeAk) 


(59) 


The  final  consideration  we  shall  give  to  Stansfield’s  two-dimensional,  radio  DF  analysis  is  directed 
toward  reproduction  of  Stansfield’s  expression  for  the  root-mean-square  error  p2  to  be  expected  in  the 
DF  position  fixes  in  which  the  statistics  are  defined  by  equation  (33)  and,  equivalently,  (47).  Given  a 
position  fix  specified  by  the  coordinates  (aq.ToK  the  absolute  positioning  error  is  defined  as 

e2  =  (V  +  yo2)'4  (60) 


The  mean  square  error  is  then 


OO  CO  CO 


P22  = 


1 


2nab 


SSS(x°2+ y°2) 


—  oo  — oo  — oo 


exp 


V 

a2  b2 


dx  dy 


(61) 


Evaluation  of  the  integral  by  expanding  the  integrand  gives 


Pi 


8(a2  +  b2) 


7 r 


OO  CO 

j  f 


u2e~"'du\ 


(62) 


Thus: 


p2  =  (a2  +  ft* 


Use  of  equations  (43)  and  (46)  with  equation  (63)  gives 

A  +  ft 
\(i  —  v2 


(63) 


Pi  = 


(64) 
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SECTION  3 

THE  THREE-DIMENSIONAL  EXTENSION  OF  STANSFIELD'S  THEORY 
3.1  THE  THREE-DIMENSIONAL  DF  ENCOUNTER 
3  I  I  INITIAL  CONCEPTS 

Figure  3-1  illustrates  the  general  distribution  of  the  participants  and  identifies  most  of  the  param¬ 
eters  used  to  characterize  the  three-dimensional  DF  encounter.  It  should  be  understood  that,  in  gen¬ 
eral.  the  DF  sites  may  be  placed  anywhere  in  the  three-dimensional  space  surrounding  the  actual  target 
position  at  O. 

As  in  the  two-dimensional  case,  it  is  the  objective  of  the  network  of  DF  sites  to  make  BL  meas¬ 
urements  and  process  these  measurements  in  an  operations  center  to  determine  a  best  estimate  of  the 
location  of  the  target  transmitter.  However,  contrary  to  the  approach  employed  in  the  two-dimensional 
case,  most  of  the  fundamental  quantities  of  interest  will  be  represented  as  vectors. 

This  could  have  been  done  in  the  two-dimensional  case,  but  it  was  desirable  to  follow  Stansfield’s 
procedure  using  scalar  quantities  only  in  order  to  facilitate  comparison  of  Section  2  with  Stansfield’s  ori¬ 
ginal  paper.  In  addition  to  considering  the  three-dimensional  theory,  this  section  will  present  the  most 
important'  of  Stansfield’s  results  as  special  cases  of  the  three-dimensional  theory  applied  to  a  planar 
space.  Where  appropriate,  these  will  be  presented  in  vector  notation. 

The  points  of  the  "encounter  volume"  throughout  which  the  DF  sites  are  distributed  are  labeled 
via  a  fixed,  right-handed  Cartesian  coordinate  system  whose  origin  is  coincident  with  the  actual  position 
of  the  target  transmitter.  It  is  assumed  that  each  DF  site  employs  some  form  of  directional  antenna 
system  from  which  a  unit  "bearing  vector"  ( B )  along  a  line  toward  the  target  can  be  determined.  The 
situation  for_the  j’th  DF  site  is  shown  in  detail  in  Figure  3-1.  The  BLjo  the  actual  transmitterjocation 
is  the  line  JO  which  is  represented  by  the  vector  D,.  The  length  of  J O  is  the  magnitude  of  D,  and  is 
represented  by  the  symbol  D,.  A  similar  convention  will  apply  for  all  vector  quantities;  e.g..  />,  is  the 
magnitude  of  the  vector  pn  thus  p,  =  p,pr 

As  in  the  two-dimensional  DF  encounter,  one  of  the  DF  sites  in  the  three-dimensional  DF 
encounter  is  considered  to  be  an  operations  center  at  which  the  BL  data  from  each  DF  site  are  gathered 
and  processed  in  order  to  determine  the  best  estimate  for  the  location  of  the  target  transmitter.  If  the 
point  5  with  coordinates  (x.v,:  )  in  Figure  3-1  is  hypothetically  offered  as  the  location  of  the  target 
transmitter  by  the  personnel  of  the  DF  operations  center,  then  the  perceived  error  associated  with  the 
target  transmitter  bearing  measurement  B,  from  the  j'th  DF  site  is  the  vector  q,.  The  actual  bearing 
error  in  the  j’th  DF  site’s  measurement,  which  is  unknown  to  the  personnel  of  the  DF  operations 
center,  is  the  vector  pr  The  vector  r  is  a  position  vector  describing  the  point  S  relative  to  the  fixed 
coordinate  system  with  its  origin  at  (>,  i.e.,  the  components  of  T  are  (x.v.: ).  Since  it  is  not  necessarily 
clear  from  inspection  of  Figure  3-1..  it  will  be  emphasized  here  that  the  vectors  Br  D..  pr  qn  and  ~r 
need  not,  in  general,  be  coplanar. 

3.1.2  THE  CONSTRUCTION  OF  p,  AND  q. 

We  embark  here  upon  a  digression  required  in  order  to  lend  credence  to  concepts  employed  later 
in  the  fundamental  statistical  characterization  of  the  three-dimensional  DF  encounter.  As  indicated  in 
Figures  3-1  and  3-2,  we  have  thusjar  assumed  the  lines  OP  and  S(2to  be  constructed  so  that  they  are 
perpendicular  to  the  line  along  JQ,  i.e.,  the  vectors  p,  and  q,  are  orthogonal  to  the  vector  Br  This 
approach  is  consistent  with  the  construction  of  analogous  quantities  in  Stansfield’s  two-dimensional  DF 
theory  and  is  sensible  when  attempting  a  logical  extension  of  that  theory  to  three  dimensions.  How¬ 
ever,  the  criterion  of  constructing  p,  and  q,  orthogonal  to  the  vector  B ,  is  otherwise  arbitrary.  Indeed, 
it  will  soon  be  convenient  to  think  of  /">,  as  being  orthogonal  to  D,  and  q ,  as  being  orthogonal  to  D,  +  r 
as  illustrated  in  Figure  3-3.  It  is  the  objective  of  this  subsection  to  demonstrate  that  the  assumption 
that  the  linear  dimensions  of  the  volume  of  unceitainty  associated  with  the  target  transmitter’s  location 
are  small  when  compared  with  the  distances  D,  (as  per  Stansfield’s  assumption  in  two-dimensions) 
leads  to  the  same  essential  analytical  conclusions  whether  we  assume  p ,  and  q,  orthogonal  to  B,  or  that 
p,  is  orthogonal  to  D.  and  q,  is  orthogonal  to  D ,  +  7; 


|  Figure  3-3  A  Magnified  View  of  the  Coordinate  Origin  to  Determine  q , 

|  When  Tjj  is  Orthogonal  to  D,+r  and  Pj  is  Orthogonal  to  £>, 
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As  in  Stansfield's  two-dimensional  DF  theory,  we  will  ultimately  employ  expressions  for  each  perceived 
bearing  error  q ,  in  terms  of  its  associated  actual  bearing  error  p ,  and  the  components  of  the  vector  r. 
We  shall  find  the  appropriate  expression  by  first  expressing  q,  in  terms  of  p,  and  r. 

First,  consider  Figure  3-2  and  the  case  wherein  p ,  and  q,  are  orthogonal  to  Br  From  the  figure 
one  can  sec  the  relation 


i 

i 

\ 


1 

I 

I 


so  that  solving  for  q,  one  has 


( q ,  +  r)  —  p  =  (q,  +  r)  •  B,\B, 


q,  =  P,  -  r  +  (r  ■  B,)  B, 


(65) 


(66) 


Now,  as  per  Stansfield’s  as.  umption,  we  shall  stipulate  that  the  BL  errors  will  be  sufficiently  small  so 
that  we  may  substitute  D,  for  5,  in  equation  (66)  without  significantly  altering  qr  and  thus  we  write 


?,=  /»,-  r  +  (r  •  b,)D, 

From  this  expression  we  find,  upon  taking  the  scalar  product, 

q,2  =  q,  ■  <7,  =  (?  -  p,)  ■  (?  -  pt)  -  •  by 


(67) 


(68) 
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Now,  let  us  attempt  the  same  thing  except  let  us  begin  by  considering  Figure  3-3  and  the  case  where  p, 
is  orthogonal  to  D,  and  cjj  is  orthogonal  to  D,  +  r.  From  the  figure,  we  see  that 

(7  +  q)  -  p,  =  a  B,  (69) 

Taking  the  scalar  product  with  B,  on  both  sides  gives 

a  =  j(?  +  q)  -  p,j  ■  B,  (70) 

so  equation  (69)  can  be  written  in  the  symmetrical  form 

Qj  ~  (q,  ■  Bj) B,  =  (p,  -?)-  [(p,  -  7)  •  £,] B,  (71 ) 

Equation  (71)  has  a  solution  of  the  form 

qj  =  Pj  -  r,  +  k,B,  (72) 

where  k,  is  an  arbitrary  constant.  However,  we  also  have  the  condition  that  q ,  is  orthogonal  to  D,  +  7 
expressed  as 

( D ,  +  7)  ■  q,  =  0  (73) 

Using  equation  (72)  in  equation  (73)  leads  to  an  expression  for  k,  given  by 

h  _7  ■(b,  +  7-pi) 


(Dj  +  7)  •  Bj 


Thus,  the  exact  solution  for  q,  is  given  by 


7  •(£>,  +  r  -  p  )  j. 

q.  p—r+  - ■=-!■ - —  B. 

( D,  +  7)  B ,  J 


Comparison  of  equation  (75)  with  equation  (67)  gives  the  impression  that  the  two  are  distinct 
which  is  the  case  in  general.  However,  we  may  rewrite  equation  (74)  as 


D,  + 

r  aa 

— -  r 

■  B, 

f 

D< 

1 

r  ~  P, 

r 


4  }  b> 


so  that  when  an^  (D^B,),  we  see  that 

ki~r  Dj 


so  that  equation  (75)  becomes 
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‘I,  =  ft  -  '•  + 


7  •  />,]/>, 


which  is  identical  to  equation  (67). 


What  we  have  demonstrated  here  is  that  when  the  linear  dimensions  of  the  volume  of  uncertainty 
of  the  transmitter  position  fix  are  small  relative  to  the  distances  £>,,  construction  of  the  vectors  p ,  and 
q,  orthogonal  to  the  vectors  D,  and  D,  +  7,  respectively,  yields  the  same  approximate  analytical  rela¬ 
tionship  between  q,  and  p,  as  when  ft  and  q,  are  constructed  orthogonal  to  Br 

Equation  (68)  plays  the  same  important  role  in  the  three-dimensional  DF  theory  as  that  played  by 
equation  (8)  in  Stansfield’s  two-dimensional  theory.  Although  equation  (78)  is  written  as  a  function  of 
the  vectors  ft,  7,  and  DJ ,  it  can  also  be  expressed  in  terms  of  BJf  Dr  and  7,  since  p ,  is  completely 
specified  if  fi,  and  b,  are  known.  Indeed,  when  p,  is  constructed  orthogonal  to  BJX  we  have  the  rela¬ 
tions 

p,  ■  Bj  =  0  (79) 


(tf,Xft)  -6,  =  0 


which  when  solved  yield 


(BjXp)  •  (D,xB,)  =  |[(£>,xfl,)  •  (b,xB,) ]*|-  H, 


(1  -  B,x  2)  Dn  +  ( B,xBn)DJX  +  (B,xBr)D,: 

H, 


(BnB,x)DJX-(  1  -  Bn)2D,x  +  (BlxBj;)D,: 


(B,:BJX)D,X  +  (B,:Bn)Dn  -  (1  -  Br  2)  D,. 


where 


H,  =  (Bj.  DJt  -  Bn  Djb2  +  (BIXDJ:  -  B,:DJXV  +  (Bn  Dj:  -  B,:DJV)' 


and  (ft, -ft, .P,:)  are  the  components  of  p ,  while  ( BIX,BIV,BJ: )  and  (/)„,£>,,  ,Dr)  are  the  components  of 
Bj  and  D,,  respectively.  As  before,  the  vertical  bars  around  terms  in  equations  (81)_and  (85)  indicate 
that  the  absolute  value  of  each  term  is  to  be  used.  Thus  p,  is  known  when  Bj  and  D,  are  known  and, 
since  the  magnitude  of  p,  is  given  by 


then  the  components  of  ft  =  ft  ft  are 


I .  ...y  .  - .  • 
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( Pj ),  =  -(1  -  B,r)  DnD ,  +  (BnBp)DnD,  +  (B/x  Bi:)  Di:Dj  (87) 

(p,),  =  (Bp  B)X)DrD,  -  (1  -  B^DjyDj  +  (BJVBJ:) DrDt  (88) 

(Pj):  =  (Bj.BjJ  D,x  D,  +  (B^BJD.D,  -  (1  -Bj:2)Dj:Dj  (89) 

Now,  when  Pj  is  constructed  orthogonal  to  D ,,  we  have  the  relations 

p,  ■  D,  =  0  (90) 

<p,xZ),)  fl;  =  0  (91) 

(D,xA>  '  (bjXB,)  =  |  [(A><4)  •  <AX  A>)|  =  H,  (92) 


which  when  solved  yield  the  components  of  p,  as 


Pj,  = 


(1  -  Dn2)B,,  -  (DnDlt.)Bjy  -  (DJXD,:)Bi: 


H, 


Pn  = 


-(DjxD„  )Bjk  +  (1  -  ZV)fl„  -  (D,,DJ:)B,: 


H, 


Pj:  = 


~  (A.-AM,  +  n  - 


(93) 

(94) 

(95) 


where,  as  before,  //,  is  given  by  equation  (85),  The  components  of  p,  are  therefore  given  by  [since 
equation  (86)  is  still  applicable] 

(96) 

(97) 


(Pj)x  =  (1  -  2)  B,x  Dj  -  (C,,Dp)BltD,  -  (DjkDj:)  Bi:Dj 

(p;),  =  -(D,XDIV)  B,x  D,  +  (1  -  Dp2)BJXD,  -  (DnDj:)B,:D, 
(p,):  =  -(DJ:D,X)BJXD,  -  (D,:DJt)Bn  D,  +  (1-Z>,;2)  BrD, 


(98) 


fxrrFWi: v&rszxURF'-  T15SB*  W'ST* JTf  *r*^'  /'^j’SrV^ 
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3.  1.3  THE  FUNDAMENTAL  STA  TISTIC'AL  CHARACTERIZA  TION  OF  THE  THREE-DIMENSIONAL 
DE  THEOR Y 

In  Stansfield's  two-dimensional  DF  theory,  the  actual  bearing  errors  p,  were  assumed  to  be  the 
values  of  random  variables  R,  that  were  normally  distributed  with  zero  mean.  In  the  three-dimensional 
case,  our  bearing  errors  are  represented  by  the  vectors  p,  which  vary  in  some  stochastic  fashion  from 
one  DF  location  attempt  to  another.  In  contrast  to  the  two-dimensional  theory  wherein  only  azimuthal 
bearing  errors  were  considered,  we  must  consider  both  azimuthal  and  elevation  bearing  error  com¬ 
ponents  in  the  three-dimensional  theory.  In  order  to  facilitate  the  characterization  of  this  slightly  more 
complicated  circumstance,  we  imagine  the  construction  of  a  right-handed  triad  of  orthogonal  unit  vec¬ 
tors  at  the  point  Dfor  each  DF  site.  For  each  DF  site,  these  unit  vectors  are  labeled  «,  j3,  and  y;  and 
the  specific  triad  for  the  j'th  DF  site  is  thus  composed  of  the  vectors  aJt  pr,  y ,.  These  vectors  define  a 
coordinate  system  (<x,p,y) ,  in  terms  of  which  the  actual  bearing  error  p ,  will  be  resolved  for  analytic 
convenience. 

Figure  3-4  is  representation  of  the  {a. Pry)  coordinate  axes  associated  with  the  j'th  DF  site.  The 
vector  «,  is  collinear  with  the  vector  D,  or  (£>,)  but  opposite  in  direction.  The  unit  vectors  p,  and  y, 
are  constructed  so  they  lie  in  a  plane  perpendicular  to  the  vector  and  such  that  p ,  lies  in  the  .v-y 
plane  while  the  projection  of  y ,  onto  it.  is  always  greater  than  or  equal  to  zero.  Mathematically,  these 
conditions  may  be  expressed  by  the  equations 

a ,  ■  y j  =  0  (99) 

(a,xy,)  •  f/_.  =  0  (100) 

(a,Xy,)  ■  (ai'XU;)  =  [«2V  +  a2]  (101) 

Solving  for  y,  from  these  in  terms  of  a,  and  using  P,  -  (y,  x  «,)  and  a,  =  -  kj ,  we  have 

a,  =  ~(D,x)ux  -  (Dn)ut  -  ( Di:) it.  (102) 


-  =  -(Plv)iL  + 

(D,x2  + D„2V 


-( D,x  Dj:)  f/,  -  ( Dn  DI:)  m,  +  ( Dn  2+ Dn  2)  uz 


<D,2  + D,2)* 


These  expressions  will  be  of  value  in  later  discussions. 


Figure  3-5  shows  how  the  actual  bearing  error  vector  p,  can  be  resolved  in  the  ( a,p,y ),  coordi¬ 
nate  system.  Since  we  shall  choose  to  construct  p,  orthogonal  to  the  vector  Dr  Pj  lies  in  the  plane  of 
y j  and  pr  A  fundamental  assumption  of  the  three-dimensional  analysis  pursued  here  is  that  the  eleva¬ 
tion  and  azimuth  components  of  (/),)y,  and  (/7,)p,  respectively,  are  the  values  of  normally  distrib¬ 
uted  random  variables  with  zero  means  and  equal  variance;  i.e.,  (P,)y  has  the  probability  density 


P  [(/>,),]  = 


.  : W 

‘P  2  trj2 
(TjUttY2- 


and  ( PA  a  has  the  probability  density 
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P 


exp 


-</>,) 2 /. 

2<0? 


<r  jO.it)'' 


(106) 


where 


(Tj  =  (<7,)yZ),  —  (<r,)pD(  =  w;D,  (107) 

Thus,  the  quantity  p,  which  is  the  magnitude  of  the  vector  />,  and  is  given  by  the  expression 

P,  =  [(£,)y2  +  ( P/>/3 2]  (108) 

is  the  value  of  a  random  variable  P,  that  is  Rayleigh  distributed  (References  7  and  8  );  i.e.,  P,  has  the 
probability  density 


pip ,1  = 


2(/>,)exp 

2tr 


(109) 


Also,  the  angle  4>,  is  the  value  of  a  random  variable  ‘F,  that  is  uniformly  distributed;  i.e.,  4),,  has  the 
probability  density 


/*(<£,]  = 


_L 

2  TT 


(110) 


These  results  characterize  the  actual  bearing  error  vector  p,  =  ptp,  as  being  composed  of  two  dis¬ 
tinct  stochastic  factors,  a  magnitude  p,  that  is  Rayleigh  distributed,  and  a  unit  vector  /;,  which  assumes 
random  orientations  about  the  line  along  D.  with  uniform  likelihood.  We  now  have  a  statistical 
representation  of  our  actual  bearing  emr  p ,  analogous  to  the  statistical  characterization  given  to 
Stansfield’s  actual  bearing  error  by  equation  (3). 

Ai  this  point  it  is  possible  to  characterize  the  set  of  perceived  bearing  errors  {</,}  resulting  from  a 
hypothesized  target  transmitter  location  S  as  a  statistical  sample  whose  likelihood  of  being  witnessed  is 
maximum  for  the  best  estimate  of  the  target  transmitter's  location. 

This  is  a  procedure  analogous  to  Stansfield's  approach  for  finding  the  best  estimate  of  the  target 
transmitter’s  location.  However,  we  will  circumvent  *his  lengthier  alternative  and,  instead,  invoke  the 
least  squares  criterion  directly.  Thus,  we  shall  stipulate  that  the  personnel  of  the  DF  operations  center 
select  as  their  DF  fix  the  point  which  minimizes  the  sum  of  the  squares  of  the  relative  errors,  where 
the  relative  error  for  the  j'th  DF  site  is  given  by 


Q, 

'  2  it,D, 


(111) 


analogous  to  the  definition  of  relative  error  giver,  for  the  two-dimensional  case.  The  2  in  the  denomi¬ 
nator  appears  because  Q ,  will  be  Rayleigh  distributed  in  the  three-dimensional  case  as  opposed  to  nor¬ 
mally  distributed  in  the  two-dimensional  case.  The  sum  of  the  squares  of  the  relative  errors  now  has 
the  form 
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ZIRijJjL 

2Ur2D2) 


analogous  to  that  of  equation  (6). 

Now,  with  the  aid  equation  (68),  equation  (112)  can  be  written  as 


(r  -  p ,)  •  (/  -  p,)  -  (r  •  D,)2j 


2  (a.D,)2 


^■a,  [ p,2  +  (x2  +  /  +  r)-2(.\u,  +  vpn  +  zpr)  -  (xDn  +  yD,x  +  r/),.)2] 

>  i - - - ^ - 1  (114) 

2  ((TjD,)2 


Equation  (114)  is  analogous  to  the  result  obtained  by  Stansfield  as  expressed  in  equation  (9),  Just  as 
equation  (9)  was  crucial  to  Stansfield’s  further  development  of  the  two-dimensional  DF  theory,  so  is 
equation  (114)  crucial  as  a  prerequisite  to  the  development  of  the  substantial  analytical  results  of  the 
three-dimensional  DF  theory.  Beginning  with  equation  (114),  we  will  derive  closed-form  expressions 
for  the  DF  fix  coordinates  ( x,y,z )  for  any  single  attempt  to  locate  the  target  transmitter,  once  the  set  of 
actual  bearing  errors  is  known  ~pp  and  we  will  describe  the  statistical  features  of  the  triple  (.v,y,z)  when 
the  actual  bearing  errors  are  characterized  as  random  variables. 

3.2  THE  THREE-DIMENSIONAL  DF  ANALYTICAL  RESULTS 

3. 2. 1  EXPRESSION  FOR  THE  DF  FIX  COORDINA  TES 

We  will  choose  as  the  best  estimate  of  the  target  transmitter’s  location  that  point  S  with  coordi¬ 
nates  (x.y,z)  that  minimizes  A  as  given  in  equation  (114).  In  order  to  find  closed  form  expressions  for 

the  coordinates  of  5,  we  evaluate  4^-,  4^-,  and  —  and  set  each  to  zero  in  order  to  produce  three 

ox  dy  dr 

equations  in  three  unknowns  and  z  The  three  resultant  equations  are 


<M>) 

(o-;D,)2 


(1  ~  Dn2) 
( (T  ,Dt )2 


-(£>„£,) 

Ur,D,)2 


~(D„D,:) 

(o-,Z>,)2 


{r,p„) 

l.r,D,)2 


(cr,D,)2 


(l  -  D,r) 

Ur,D,)2 


-(D„DJ 

( <t,D ,)’ 


(P,P,A 

(<t,D,)2 


-(*w 

((TjD,)2 


.  v  : Ml  .  ; 

(a  ,Dt)2 


\p[  (1  -  D,;2) 

4*  ( <r,D ,)2 


Next  we  define  a  set  of  fundamental  parameters  as  follows: 
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<„8, 

I 

fl  =  y[iL;A2i| 

M  Ur.D,)2 

'-X[^ 

/ 

y  -<P„o.)|  3) 

/ 

Now,  solving  equations  (115)  through  (117)  for  .v;  yand  z  gives 

<*  -  + <*  -  wl + *  -  ”'I|  W 

x  — - 1 - { - L - 1 - '■ -  (124) 

[X(/x|  -  £2)  +  v(^Tj  -  1>€)  +  if(t>£  -  T)/x)J 

+ <*  -  ”:)I[  w| + -  «>X(  w 

V-  - 2 - ? - - - i - - -  (125) 

jx(/z£  -  £2)  +  ^({t}  —  v£)  +  T)(vl,  -  T)ju)J 

("{  ~  +  -  *{,L[  c^r] +  ~  ‘'1)E(<^r 

2  —  - ! - 1 - -  (126) 

[X<M*  -  &  +  -  v&  +  >j(K  -  t,m)| 

Equations  (124)  through  (126)  are  the  coordinates  of  the  best  estimate  of  the  location  of  the  target 
transmitter  as  determined  by  the  personnel  of  the  DF  operations  center  for  a  given  set  of  actual  bearing 
errors  Ip,  =  p.pJX ux  +  p,Pjyux  +  P,P,:u .-)■  These  are  analogous  to  Stansfield’s  equations  for  the  fix  coor¬ 
dinates  given  by  equations  (10)  and  (11).  Indeed,  if  we  apply  these  equations  to  the  case  where  all  of 
the  DF  encounter  participants  lie  in  the  x-y  plane  and  we  set  the  elevation  error  to  zero,  then 
r?  =  {  =  Pi;  =  0,  and  equations  (124)  through  (126)  become 


[X/x  -  j/2] 


brX 


P, 


JUP-Pjx  ~  VP»  I 


tyidifA&jtib  &  >ki  iMi,  i 
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^,7  ^  -  -".-I 


Now,  as  per  Stansfield’s  theory,  the  vector  is  given  by  p,  =  />,/>,,  where  p,  is  the  value  of  a  normally 
distributed  random  variable  P,  of  zero  mean  (thus  it  can  assume  positive  and  negative  values  with 
equal  likelihood)  and  p,  is  a  fixed  unit  vector  very  nearly  orthogonal  to  D,;  i.e.,  (see,  for  instance.  Fig¬ 
ure  2-1) 

P,  —  —  sin  ti,ux  +  cos  (129! 

Equations  (127)  and  (128)  can  now  be  written  as 


=  1  T'  Pi  c°s  A, 

(\/A  ~  V2)  S  (( TjDj )2  ^ 


Eh  sin  Oj 
(c. t,D,)2 


=  1  T'  P.cos  Oj  _  Tp  P, sin  Oj 

y  (A/x  -  v2)  (<TjD,)2  '  (cT/D,)2 


where  vs=  —  v. 

Equations  (130)  and  (131)  are  Stansfield’s  equations  (16)  and  (17)  determined  directly  from  the 
three-dimensional  fix  equations  (124)  through  (126).  Using  equations  (93)  through  (95)  for  the  case 
when  B)X  -  Dn  =  0,  we  can  express  Stansfield’s  equations  entirely  in  terms  of  the  set  of  D,  and  the 
components  of  B,  and  D,, 


.  1  VI 

1 

( Bn  Dn  -  D,x)(p  Dl}  +  vDls) 

4 

r-* 

1 

Vj2D, 

(D,2+  D, ,2) 

1  Tl 

1 

(BnD,s  -  BnDn)(vDn  +  \DIX) 

( A/x  —  v2) 

(Dlx2+Dlt2) 

which  can  be  expressed  in  vector  notation  as 


■^2(^7  ■ b) +  -'A  ■  A>] 

/  ’■  '  ' 
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J.  2. 2  DF  FIX  COORD!  N A  TE  ST  A  TISTIC'S 

Let  us  define  the  variables  i1(  /'2,  and  by  the  equations 


r  i  = 


i->  = 


r. !  = 


(136) 

(137) 

(138) 


We  have  changed  the  subscripts  x  to  1,  y  to  2,  and  .v  to  3  for  computational  and  notational  conveni¬ 
ence;  i.e.,  plx  =  p, |,  ph  =  p,i,  p,:  =  p, 3. 


Since  p,,  Pj 1 1  Pj2,  and  p, 3  are  all  the  values  of  random  variables,  where  P,  is  Rayleigh  distributed 
as  per  equation  (109)  and  p,\,  /j)2,  and  p,2  are  all  distributed  alike  (we  will  see  exactly  how  they  are  dis¬ 
tributed  later),  then  the  quantities  /].  r2,  and  r2  are  the  values  of  random  variables  we  will  denote  as 
R |,  R2,  and  R2.  Note  that  the  random  variables  V\  =  PjP,\ ,  l'2  =  P,Pl2x  and  Kj  =  Pj Pl2  are  all  simi¬ 
larly  distributed.  Inspection  of  equations  (124)  through  (126)  reveals  that  the  fix  coordinates  .v,  .v,  and 
c,  are  linear  functions  of  r, ,  r2,  and  r}  and  thus  the  values  of  random  variables  we  will  denote  as  V.  Y, 
and  Z  It  is  our  objective  in  this  section  to  determine  the  joint  probability  density  for  the  random  vari¬ 
ables  -V,  K  and  Z  whose  values  are  the  fix  coordinates  for  any  given  attempt  to  locate  the  target 
transmitter. 


We  will  begin  by  assuming  that  the  number  of  DF  sites  n  is  large  enough  that,  by  virtue  of  the 
central  limit  theorem  (Reference  5).  the  random  variables  R\,  R2,  and  R2  are  normally  distributed. 
The  joint  probability  density  for  R\t,  R2 ,  and  R2  is  then  given  by  (References  6  and  9) 


p('i.'2.'j) 


exp 


y([/?n \r'‘)[R} 


(277)3/2(Det[C))'- 


(139) 


where  (Cl  is  the  covariance  matrix,  [C  ’]  the  inverse  of  the  covariance  matrix,  [R]  is  the  deviation 
matrix,  and  [/?']  is  the  transpose  of  the  matrix  (/?].  The  determinant  of  the  matrix  [C]  is  written  as 
Det  [Cl. 

The  elements  of  (Cl  are  given  by 


C,„„  =  Cov  (R„.RJ 


(140) 


where  n  =  1,  2,  or  3  and  m  =  1,  2,  or  3;  i.t.,  C.„„  is  in  the  mth  row  and  nth  column  of  the  3-by-3 
matrix  (Cl.  The  matrix  [/?!  is  a  column  vector  given  by 


IR] 


1 1 
>2 


~(K  1} 

- 


(141) 


where  [R\),  (R2),  and  (Ri)  are  the  mean  values  of  the  random  variables  Rt,  R2 ,  and  R'„  .-.-j;  ■  mvc-ly. 
To  determine  the  explicit  form  of  the  joint  density  of  equation  (139),  one  must  find  [R\).  and 
(R  ,)  and  the  covariance  elements  C„„„  n  =  1,  2,  3  and  m  =  1,  2,  3. 

From  the  definitions  of  the  mean  and  covariance 
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/  R  \ 

\A'V 


II 


/«/ 


Ut,D,)2 


//= 1,2.3 


(142) 


and 


( nit 


[o-^aADJ2 


/j  and  tw=l  .2.3 


(143) 


It  now  becomes  apparent  that  we  must  know  the  probability  density  for  each  of  the  random  variables  in 
equations  (142)  and  (143).  We  already  know  that  Pj  is  Rayleigh  distributed  as  per  equation  (109),  and 
we  will  now  determine  the  probability  density  for  each  random  variable  P,„,  n  =  1,  2,  3.  whose  values 
are  the  components  of  the  unit  vector  pr  Before  proceeding,  however,  we  shall  impose  the  assumption 
that  P,  and  Pm,  n  =  1,  2,  3,  are  statistically  independent  quantities;  i.e.,  we  suppose  that  the  magnitude 
of  p,  and  its  direction  are  independent.  We  further  suppose  that  the  measurements  at  one  DF  site  are 
independent  of  those  at  another  DF  site.  Equations  (142)  and  (143)  can  now  be  written  as 


I  r  \  = 

\K»I 


n  and  w=l ,2,3 


(144) 


C  — 

v  nni 


/  p  i\/  p 
\ri  /\‘ i 


Plin)-{P)\P^P, 


(cr  ,D,)4 


,n  and  m—  1 ,2.3 


(145) 


Now,  notice  from  Figure  3-5  that  p ,  can  be  written  as 

f>,  =  cos  +  sin  4>jy,  (146) 

and,  since  (3,  and  y,  are  given  by  equations  (103)  and  (104)  respectively,  equation  (146)  can  be  written 
as 

P,  =  ~2  ■  |(D,j  cos  rin 

(D,\l  +  Di2-)'  I 

-  (£), i  cos  <{>,  +  Dti  sin  ,)  F/>  +  ( D,\ 2  +  Dti2)  sin  (147) 

from  which  we  see 


p,\  = 

( Da  cos  (f>,  -  D, |  D,\  sin  $,) 

(148) 

(Dr  +  Da2) 

Pa  = 

-  (D,,  cos  <$>,  +  DaD,y  sin  <t>,) 

(149) 

(D,\2  +  Da2)" 

(D,r  +  Da2)'"  sin  6, 

(150) 

Each  of  the  components  of  Pj  varies  in  a  random  manner  because  the  angle  <f>,  is  the  value  of  a 
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uniformly  distributed  random  variable  <P(  [see  equation  (110)].  The  statistical  nature  of  the  random 
variables  Pg,.  P, 2,  and  P, 3  is  more  readily  recognized  if  equations  (148)  and  (149)  are  written  as 


I’ll 


D,y  +  W 

V  +  D,y 


sin  [<)>,  +  S|] 


(151) 


P,2  = 


V  +  W 


0,1  2  +  D,22 


sin  [<^,+82] 


(152) 


where 


8t  =  tan  1 


82  =  tan  1 


(-A,^ 


Da 
Dg  Dg 


Da 


(153) 


(154) 


Thus,  the  random  variables  P„„  n  =  1,  2,  3,  have  probability  densities  of  the  form  (References  7  and 
8) 


pIpJ 


1 


in*  /  r  2  *>\ •/. 

^  ^ in  Pm  ' 


(155) 


where  F„,  is  the  coefficient  of  the  sine  factor  of  p„,  in  the  forms  of  equations  (150)  through  (152);  e.g.. 


for  Fj\ 

(/*»};  i  e., 


Dg2  +  Dg2,  D,y 


D„2  +  DJ 


.  Thus,  \pw\  ^  Fitl  .  With  equation  (155)  we  are  able  to  evaluate 


{P,n)  =  P„APJ  ‘h>„,  =  JflZ 


dRn —  =  o 
pJY- 


(156) 


a  result  that  follows  from  the  fact  that  the  integrand  in  the  second  integral  is  an  odd  function  of 
Equations  (144)  and  (145)  now  reduce  to 


{*„}  =  (> 


Cm 


From  equation  (109)  it  follows  that 


X{Pr){p„A„) 


<F,2)  =  2(cr,D,  2 


thus  we  are  left  with  evaluating  {P,nPim).  In  order  to  do  this,  we  first  express  p„,  as 

Pm  =  am  cos  <t>,  +  bm  sin  <t>, 
a  simplified  representation  of  equations  (148)  through  (150),  where 


(157) 

(158) 

(159) 

(160) 
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aJ  1  = 

(0,i 

2  + 

0,2  2)" 

-0/i 

ai2  ~ 

(Ai 

2  + 

0,22)% 

= 

0 

- 

-0,i 

0/2 

b,\  ~ 

(D„ 

2  + 

0,2  2)* 

u 

- 

-0,? 

0/3 

bji  - 

(O/i 

2  + 

0,22)'/; 

Z,/3  =  (D/12  +  D/22>* 

In  general,  then, 

PjnP/m  =  $/  ^m^nn  Qjmbjn)  COS  4>  ,  Sin  4* ,  +  bjnbjm  Sin  <j)  j 

thus  use  of  equation  (110)  gives 


lit  j  *7T 

( Pjn  Pjmj  ~  J P/n(<i> /)  PmS^  )  ~  *  /  Pjn  ($  Pjm  (4>j)  d<bj 


Therefore,  use  of  equation  (167)  in  equation  (168)  yields 

_  f  27T  2tt 


\pjnp,m'j  =  4~-  amawS  c°s2  4>,  (I <t> ,+((!,„ bl,„+aimblll)  j  cos  <t>,  sin  <b,  cl<t>j+binblinj  sin 2  <t>j  d<l>,  (169) 

i'  i\  n 


P  p  \  = 

1  ///*  />»;/ 


M/ll  U/ni  "1"  bm  b/nl 


Now,  substitution  of  equations  (159)  and  (170)  into  equation  (158)  gives 


r  = 

^  IWI 


*t~  b^bjin) 

(<r  ;Z),)2 


Explicit  substitution  of  the  a's  and  b's  from  equations  (161)  through  (166)  gives  for  the  covariance  ele¬ 
ments  (upon  changing  back  to  the  subscripts  .v.r.r  for  1,  2,  3,  respectively) 


r  =T[lLl_Dii 1 
"  M  (.A)2 


a  -  v> 

(cr,Dj)2 
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c  _Y  [Hi All 

La  Ur,!),)1 

r  =r  r 

",  c"  La  ((TiO,)2 


Q  =  Q  = 


-<A'A:> 

Ur,D,)2 

-(A,  A,) 


C  -C  - Y 


Comparison  of  the  covariance  matrix  elements  with  equations  (118)  through  (123)  shows  that  the 
coefficients  of  the  fix  coordinates  (.v, v,c  )  in  equations  (115)  through  (117)  are,  in  fact,  these  same 
covariance  elements.  Also,  equation  (157)  means  that  the  deviation  matrix  of  equation  (151)  reduces 
to 


Thus,  if  we  define  the  column  vector  [X]  as 


then  equations  (115)  through  (117)  can  be  expressed  by  the  matrix  equation 

[/?]  =  [C][T]  (180: 

This  important  result  provides  us  with  the  capability  to  use  equation  (139)  for  the  joint  probability  den¬ 
sity  of  the  random  variables  /?,,  /?, ,  and  R:  in  order  to  derive  the  joint  probability  density  for  the  DF 
fix  random  variables  X  Y,  Z.  Indeed,  since  i\,  r,,  and  r.  are  each  functions  of  the  fix  coordinates 
(x,v,z)  as  expressed  b>  (180),  then  the  joint  probability  density  for  the  fix  coordinates  is  given  by 

p(x,v,z)  =  !  p  ( r,  [  ,vj>,  c  ],/,!  .v..r, :].  r.  [  .vj .  c  ])  H  7  (/■,/•,/•-;  xf.  r )  1  (181! 

where  J  (rx,rx,r.\  x,y,z)  is  the  Jacobian  determinant  for  the  transformation  from  the  coordinates 
(rN,r,,/--)  to  the  coordinates  (.vju).  Use  of  equations  (1 15)  through  (117)  yields 


J{rx  rv  r.  x,y,:)  =  DettC] 


From  equations  (139)  and  (180)  we  have 


p  rx[x,y,z) ,ry{x,y,z) ,r.[x,y,z)  = 


expj-junicntc^ncim 

(27r)3/2(Det  [C])v‘ 


or,  since  [C^'KOlAl  =  l-T]  and  the  symmetry  of  [C]  gives  (Cr)  =  1C],  then 


-  29  - 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL,  Maryland 


exp 


-ju'Hcim 


(27r)w(Det[C))'/; 


(184) 


Use  of  the  parameter  definitions  in  equations  (118)  through  (123)  and  substitution  of  equations  (182) 
and  (184)  into  equation  (181)  gives  the  three-dimensional  joint  probability  density  for  the  DF  fix  coor¬ 
dinates  as 


p(x,y,z )  = 


A(/n£  -  £2)  +  ^({TJ  -  f£)  +  T}(l/£  -  7}/Lt) 


(2n) 


.1/2 


exp 


yfA.v2  +  py2  +  ^z2  +  2i >xy  +  2  t)xz  +  2£yz) 


(185) 


After  the  expressions  given  for  the  fix  coordinates  in  equations  (124)  through  (126),  equation 
(185)  is  the  second  significant  analytical  result  of  the  three-dimensional  DF  theory.  When  multiplied 
by  the  differential  volume  of  the  neighborhood  about  a  point  (.v.y’,:),  the  joint  density  of  equation 
(185)  gives  the  probability  that  the  DF  fix  resulting  from  any  single  attempt  to  locate  the  target 
transmitter  will  fall  within  that  neighborhood  of  the  point  (jr.y'.z).  The  quadratic  nature  of  the  argu¬ 
ment  of  the  exponential  and  the  fact  that  A,  p,  and  £  are  all  greater  than  or  equal  to  zero  indicates  that 
the  surfaces  of  constant  likelihood  are  ellipsoids  centered  on  the  actual  target  position  (the  point  O  in 
Figure  3-1). 

Equation  (185)  is,  of  course,  analogous  to  that  of  Stansfield's  two-dimensional  theory  as  given  by 
equation  (33).  Indeed,  if  one  applies  the  statistical  analysis  procedure  employed  to  generate  equation 
(185)  to  the  case  where  all  the  participants  of  the  DF  encounter  lie  in  the  x-y  plane  and  the  elevation 
error  p,:  =  0  for  each  DF  site,  then  the  matrix  [/?]  becomes  reduced  to 


(/?]  = 


(186) 


where  r,  and  r,  are  given  by  equations  (136)  and  (137).  One  also  considers  that  the  entire  random 
variation  of  the  bearing  errors  is  embodied  in  normally  varying  random  magnitudes  Pr  the  unit  vector 
components  p„,  being  deterministically  related  to  the  vectors  Dr  The  random  variables  Rx  and  /?,  are 
still  normally  varying,  so  the  joint  probability  density  for  the  random  variables  7?,  and  /?,  is  given  by  a 
modification  of  equation  (139)  for  the  two-dimensional  case;  i.e.. 


/  \  _  exp[  —  ((/? r] (C  *][/?))] 

P  xJy  ~  (27r)'/;(Det  lC])'/: 


(187) 


The  covariance  matrix  is  now  a  2-by-2  matrix  whose  elements  are  given  by  equations  (172),  (173),  and 
(175),  where  it  must  be  understood  that  the  vector  D ,  has  only  the  components  and  in  this 
two-dimensional  case.  We  see  then  that  if  we  now  define 

[X]  =  [^]  (188) 

then  we  have  an  analogous  result  in  two-dimensions  to  that  given  by  equation  (180)  for  three- 
dimensional;  i.e.. 


[R]  =  [C][,V] 
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Substitution  of  equation  (189)  into  equation  (187),  use  of  the  parameter  definitions  from  equations 
(118),  (119),  and  (121),  and  multiplication  of  the  result  by  the  Jacobian  determinant  for  the  variable 
change  from  the  (/\.r,)  space  to  the  (,v,.r)  space  gives 


P  (•'.)') 


—  yfA.v2  +  2  vxy  +  /jlv') 


(190) 


This  is,  of  course,  Stansfield’s  result  and  is  identical  to  equation  (33)  [recall  that  u  as  defined  by  equa¬ 
tion  (121)  is  opposite  in  sign  to  Stansfield’s  v  as  defined  by  equation  (14)], 


3.2.3  FURTHER  CONS/DERA  TIONS  OF  THE  JOINT  DENSITY  OF  THE  DF  FIX  COORDINA  TES 

Now  that  it  has  been  shown  that  the  surfaces  of  constant  likelihood  for  the  DF  fix  of  a  three- 
dimensional  DF  encounter  are  ellipsoids  centered  on  the  actual  target  transmitter  location,  it  is  useful 
to  deuve  ancillary  results  which  allow  us  to  describe  the  dimensions  and  orientation  of  these  ellipsoids 
relative  to  our  reference  coordinate  system  as  well  as  to  compute  the  probability  that  a  DF  fix  will  fall 
within  the  region  bounded  by  these  surfaces. 

From  equations  (181),  (182),  and  (184)  we  have 


p(.v,v,r) 


[Det  (Cl}*expl-jU7)  l  C]{X)) 
(2n)in 


(191) 


which  is  an  alternate  form  of  equation  (185)  for  the  joint  probability  density  of  the  DF  fix  coordinates. 
The  scalar  quantity 

£>=  UMlClM  (192) 

which  apart  from  a  factor  of  (-  Vi)  is  the  argument  of  the  exponential  in  equation  (191),  is  a  quadratic 
form  because  of  the  fact  that  [C]  is  a  symmetric  matrix.  Since  [C]  is  also  real,  it  can  be  diagonalized 
(Reference  [10])  by  an  appropriate  transformation.  In  particular,  there  is  an  orthogonal  transformation 
matrix  [A]  which  transforms  the  coordinates  (aj’.c)  of  a  point  F  in  the  encounter  space  relative  to  our 
reference  coordinate  system  to  the  coordinates  (awp.-o)  of  the  same  point  ^relative  to  a  coordinate 
frame  whose  origin  is  stili  at  the  actual  target  position  but  whose  axes  are,  in  general,  rotated  relative  to 
the  axes  of  the  reference  system.  Thus,  there  is  a  matrix  [A]  such  that 

U0!  =  \A][X]  (193) 

where 

[Xn]  =  v0  (194) 

-o 

Since  the  transpose  of  an  orthogonal  transformation  matrix  [.4]  is  also  its  inverse,  then  equation  (193) 
implies  the  results 

\Xr]  =  [(T0),][.4).[A']  =  U'H.y  (195) 

so  equation  (192)  becomes 


q  =  [<*„)  nunciunuo] 


(196) 
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where  we  require 


l/lllC'JU']  =  ID] 


where  [D]  is  a  diagonal  matrix.  If  we  denote  the  matrix  element  in  the  j'th  row  and  k'th  column  as 
[Dl,*,  then  equation  (191)  becomes 


pU.ro. :o)  =  {DettCllV(2^3^^— j  exp  {- ^-(ao2[ O] , ,  +  voMZ)l22+r02[Z>l33)}  (198) 


J(x,y,z\x0  y0  :a)  -  1  (199! 

for  orthogonal  transformations  to  rotated  coordinate  systems  (with  no  inversion),  then  we  may  write 
equation  (198)  as 


,  ,  _  {Pet  [Cl}/i 

PU o.Jo.V  ~  (2tt)3/2 


1  a-02  y02  z02 

exP  -y  -y- + -3- + -r 


where  to  allow  comparison  with  Stansfield’s  results  we  have  defined 


=  [Dm 


-j  =  lD]j,  (203) 

Using  the  definitions  in  equations  (201)  through  (203)  and  taking  the  determinant  of  both  sides  of 
equation  (197)  gives 

Det  [C]  =  -r—  (204 

a *  b  r 

so  equation  (200)  becomes 


p(xQ.y0,z0)  =  — 


1  ao2  yo 2  z02 

exTI  T+-F  +  ~ 


(.lir)i/2(abc) 


The  result  of  equation  (205)  is  analogous  to  that  of  equation  (47)  from  Stansfield’s  two-dimensional 
theory.  The  joint  probability  density  p(.v0.^o-2o)  f°r  the  random  variables  X0,Y0.Z0  determines,  when 
multiplied  by  a  differential  volume,  the  likelihood  that  the  DF  fix  coordinates  resulting  from  any  single 
attempt  to  locate  the  target  transmitte:  will  fall  within  a  differential  neighborhood  of  the  point 
UoJb.zo).  Thus,  it  is  possible  to  find  a  coordinate  system  with  the  actual  target  location  at  its  origin 
and  relative  to  which  the  surfaces  of  constant  likelihood  are  ellipsoids  expressed  in  canonical  form;  i.e., 
each  of  the  locus  of  points  for  which  the  argument  of  the  exponential  in  equation  (205)  is  a  constant 
value  k 2  satisfies  the  equation 
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(206) 


Equation  (206)  is  analogous  to  equation  (34)  discussed  in  association  with  Stansfield's  elliptic  con¬ 
tours  of  equal  likelihood.  It  was  subsequently  shown  in  equations  (44)  and  (45)  that  the  quantities  — , 

U‘ 

and  -t  could  be  defined  in  terms  of  the  fundamental  parameters  K,  /a.  and  «>.  In  order  to  accomplish 
b~ 

a  similar  sort  of  result  for  the  three-dimensional  theory,  the  following  identity  was  considered: 

Ul[C]Ur]  -  Ul/tllM7]  =  ID)  -  /[l]  (207) 

where  (1]  is  the  identity  matrix.  The  validity  of  equation  (207)  follows  from  equation  (197;  Using 
the  distributive  property  of  matrix  multiplication  over  addition,  equation  (207)  can  be  written  as 

lA]{[C]  -  r[lJ)Mrl  =  \D)  -  /111  (208) 

and  taking  determinants  on  both  sides  gives 

Det{[C']-/[l}}  =  n([Z>]„-/)  (209) 

Equation  (209)  is  the  equation  for  the  three-dimensional  theory  that  is  analogous  to  equations  (44)  and 
(45)  from  the  two-dimensional  theory.  Equation  (209)  shows  that  the  diagonal  elements  of  (£>]  are  the 
roots  of  the  third-degree  polynomial  in  r  on  the  left-hand  side  of  the  equation.  Once  the  elements 

1  1  1 

[£>].,,  j  =  1,  2,  3,  are  found,  then  equations  (2.01)  through  (203)  yield  —r,  ~r,  and  -r.  Indeed, 

a 1  b~  r 

when  set  equal  to  zero,  the  left-hand  side  of  equation  (209)  is  the  characteristic  polynomial  of  the  prob¬ 
lem  and  is  expressed  in  terms  of  the  fundamental  parameters  of  equations  (118)  through  (123)  as 


r*  +  (/x  +  \  +  £ )  r  —  (v2  +  7)2  +  —  fit;  —  fj.  k  —  Af)  r 

—  —  £2)  +  i/(£tj  —  v£)  +  —  Tj j/x)  =  0  (210) 

As  mentioned,  each  of  the  matrix  elements  [£>}„  is  set  equal  to  one  of  the  roots  of  the  characteristic 
polynomial  in  equation  (210),  one  root  is  associated  with  only  one  matrix  element  at  a  time.  However, 
when  there  are  distinct  roots,  there  are  six  unique  ways  in  which  the  matrix  elements  1  D]u  can  be 
assigned  values.  This  circumstance  merely  reflects  the  fact  that  there  are  six  unique  ways  in  which  the 
axes  of  the  reference  coordinate  system  can  be  rotated  and  aligned  with  the  principal  axes  of  an  equal 
likelihood  ellipsoid  so  as  to  yield  a  reference  system  relative  to  which  the  equation  of  the  locus  of  the 
ellipsoid  in  canonical  form.  For  most  applications,  the  ambiguity  ip.  the  assignment  of  values  to  the  ele¬ 
ments  f/7],,  will  be  of  no  consequence,  although  it  is  important  for  the  analyst  to  understand  how  his 
assignment  of  these  values  affects  the  characterization  of  his  particular  problem. 

Once  the  matrix  ( D ]  is  known,  one  can  also  specify  the  unit  vectors  that  define  the  axes  of  the 
(x0,y0,co)  coordinate  system  and  thereby  determine  the  orientation  of  a  given  family  of  ellipsoids  of 
constant  likelihood.  Starting  with  equation  (197),  we  find 

(C'lU?]  =  l.4n[/9)  (211) 

If,  for  a  given  square  matrix  [A/l,  we  let  [M])  denote  a  column  vector  formed  from  the  j’th  column  of 
l Ml,  then  from  equation  (21 1)  we  have 


X>1 


-  33  - 


V  (fJUPJWI  |W«  V&m  ^  }  3 


tHE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

UUREL,  MARYUND 


[lC}Mr]]„  =  [[C)l/<n;](  «  [Ar\„[D\u  (212) 

Now,  since 

ur!,/  =  [uiur);]( 

then  equation  (213)  becomes 

[icll/i'l;]  -lo]„[lilM'l;|, 

or 

[Cil/);  =  [fll„li)[/l;  (215) 

and  thus 

[[C]-U)]wm][^);-0  (216) 

Equation  (216)  is  a  matrix  representation  of  three  equations  for  the  three  unknown  l AT)\j,  t.4r] 2j,  and 
[AT]}j.  The  equation  will  have  a  nontrivial  solution  only  if 

Det[[C]— IZ>]„11]]=0  (217) 

Starting  again  with  (197)  one  can  show  that 

U][(C1— [Z>](,.ll)]Mr]  =  [Z>]-[Z>]„[11  (218) 

where  [ £>]  -  [Dlyy  [1]  is  a  diagonal  matrix  for  which 

[l/>]  -  (/)^(1][  =  0  (219) 

so  that  taking  determinants  on  both  sides  of  equation  (218)  yields  (217).  Thus,  non-trival  solutions  of 
equation  (216)  are  assured. 

We  now  have  a  way  of  finding  each  of  the  columns  of  the  matrix  [Ar],  and  thus  we  have  \AT  1 
and  the  orthogonal  transformation  matrix  l -4 ]  (which  is  the  transpose  of  [TrD.  Notice  from  equation 
(195)  that 


(213) 


(214) 


[X]  =  Ur)[T0] 


(220) 


A  unit  vector  along  the  .v0  axis  can  be  written  in  the  („y0  ..Vo.-o)  system  as 


(221) 


so  the  coordinates  of  this  vector  in  the  (x,.v,z  )  system  are  given  by  substitution  of  equation  (221)  into 
equation  (220).  A  similar  procedure  is  used  to  express  w,n  and  w-0,  yielding 


%  =  \Ar\i 
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I 

I 


(224) 


Thus,  the  axes  of  the  (.y0j'o.m))  coordinate  system  are  known,  and,  because  these  are  along  the  princi¬ 
pal  axes  of  the  ellipsoids  of  constant  likelihood,  we  know  the  orientation  of  these  principal  axes  relative 
to  the  axes  of  the  (.vj.r)  coordinate  system.  Indeed,  to  specify  the  orientation  of  any  given  ellipsoid  of 
constant  likelihood  we  now  need  only  to  determine  the  points  at  which  the  ellipsoid  intersects  each  of 
the  principal  axes,  and  this  is  easily  done  with  the  use  of  equation  (206)  since  each  of  the  matrix  ele¬ 
ments  [£>]„  is  associated  with  one  column  vector  ( AT\)  and  is  related  to  the  coefficients  of  equation 
(206)  via  equations  (201)  through  (203). 

We  are  now  in  a  position  to  consider  the  problem  of  determining  the  probability  of  the  event 
Uo\>o',z0')e  Vk  that  the  coordinates  of  a  DF  fix  (a0',voW)  will  fall  within  the  region  Vk  enclosed  by 
the  ellipsoid  of  constant  likelihood  whose  locus  is  given  by 


2  ■> 
a 2  b 2 


A2 


(225) 


A  family  of  ellipsoids  concentric  about  the  point  0  of  the  actual  target  position  is  given  by  the  equation 

~~  +  ^V-  +  =  A2,  0<  A<°°  (226) 

a2  b2  c2 

The  ellipsoid  of  equation  (225)  is  one  member  of  this  family.  The  probability  that  the  fix  point 
(Vj’o.-o')  lies  between  the  ellipsoidal  suffaces  defined  by,  A  and  A  +  dK  is,  using  equation  (205), 


-Vo. To  -n)  ilxod.vadzo  = 


(227) 


Thus,  the  desired  probability  can  be  evaluated  by 


/>(.v0'j’oW«  vt) 


4 

— —  I  w2exp[— w2]  du 

(wr  0 


(228) 


or 


P(x0',y0',z0'€  V.)  =  erf 


k 

2_ 

'/! 

A2 

(2)'/; 

7 T 

A  exp 

2 

(229) 


Equation  (229)  is  analogous  to  equation  (54),  derived  for  the  two-dimensional  DF  encounter.  We 
cannot  solve  equation  (229)  explicitly  for  k  in  terms  of  P{.\0',yQ',:n'e  Pk)  as  was  possible  for  equation 
(54)  in  two  dimensions.  However,  we  can  still  derive  results  in  the  three-dimensional  case  that  are 
equivalent  to  equations  (57)  and  (58)  which  give  the  principal  axes  intercepts  in  the  two-dimensional 
case.  Indeed,  if  one  plots  the  function  P(.\0',y0',:0'e  Vk)  as  given  by  equation  (229)  on  "probability 
paper"  (see  Figure  3-6),  then  one  can  graphically  invert  the  equation  to  solve  for  A  given  a  desired 
value  of  P(. Yo'.yoWe  Pk).  Thus,  if  one  wishes  to  determine  the  principal  axis  intercepts  of  an  ellipsoid 
within  which  the  DF  fix  point  has  a  probability  y  of  falling,  one  enters  the  axis  of  P(.v0',To'.Ve  at  T 
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and  traces  back  to  find  the  associated  value  of  A,  Av  This  procedure  is  illustrated  in  Figure  3-6.  The 
principal  axes  intercepts  of  the  ellipsoid  within  which  the  DF  fix  has  probability  y  of  falling  are  then 
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Figure  3-6  Graphical  Inversion  of  FOc'oJ’o.z'o.c  FA) 
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One  can,  in  fact,  reduce  computation  time  for  two-dimensional  DP  problems  by  graphical  inversion  of 
equation  (54)  instead  of  using  equations  (57)  and  (58)  directly. 

Of  course,  in  light  of  equations  (230)  through  (232),  the  volume  il;  of  the  region  1A  contained 
within  the  ellipsoid  is 


O.T  = 


tlbl'ky  ^ 


(233) 


If,  when  given  a  DF  fix  specified  by  the  coordinates  ( vod'o.-oD  we  define  the  absolute  positioning 
error  as 


(’)  =  (-Vo"  +  .r,) 2  +  r02). 


(234) 


then  the  mean  square  error  is  given  by 


Tl 

P3*.= 


<2„)kfcJ  + r°! + ;°!)exp 

—  CO  —  oo  —  oo 


1 

2 


-To 2  Jo2  V 

a2  b2  o2 


dx0dy0ciz0 


(235) 


Evaluation  of  this  integral  by  expanding  the  integrand  gives 


P.O  = 


(2  ir)i/2abc 


+ 


H\\  .7- 

exp(— *»“)</»»'  2 -'/2ir'J  H-2exp(-ie2)r/it 

CO  1 2  CO 

2ac  J'expf— 2 3/2  6*  j  u~exp(- u-)  < 

0  ]  0 


CO 

+  lab  J  exp(— 
o 


ir2exp(—  ir2)  dw 


(236) 


and  thus 


P}  =  U2  4 -  b2  +  t2)7' 


(237) 


This  is,  of  course,  the  three-dimensional  analog  of  equation  (63).  Using  equations  (201)  through  (203) 
we  see  that  equalion  (237)  can  be  written  as 


P?  = 


[D  ]„2  I  D)rJ  [D]jj2 


(238) 
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The  methods  of  this  section  can  be  applied  to  the  two-dimensional  DP  encounter  as  well.  In  the 
two-dimensional  case  the  matrix  (C]Js  a  2-by-2  matrix  whose  elements  are  given  by  equations  (172), 
(173),  and  (175),  where  the  vector  D ,  has  only  the  components  Dn  and  Dn.  The  characteristic  equa¬ 
tion  of  equation  (209)  then  becomes  a  second-degree  polynomial  given  by 

r2  -  (X  +  fi)  r  +  ~  v2)  =  0  (239) 


from  which  we  derive 
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{£>](!  =  -y  =  (A  +  fi)  ±  |(A  —  n )2  +  4t/21 
a  2  1  1 


[^22  =  -p-  =  j|(A+/ll)  ±  [(A  -  /J.)2  +  4*>2j 


Equations  (240)  and  (241)  reproduce  the  results  of  equations  (44)  and  (45)  as  derived  by  Stansfield. 
The  upper  signs  in  equations  (240)  and  (241)  set  the  unit  vector  wV|)  along  the  semiminor  axes  and  the 
unit  vector  w,  along  the  semimajor  axes  of  the  given  ellipse  of  constant  likelihood.  The  lower  signs 
place  f/Vi)  along  the  semimajor  axes  and  uVg  along  the  semiminor  axes  of  the  ellipse  of  constant  likeli¬ 
hood.  Indeed,  using  equations  (240)  and  (241)  to  solve  equation  (216),  and  then  using  equations 
(222)  through  (224)  gives  the  unit  vectors  of  the  axes  of  the  rotated  reference  system  (relative  to 
which  the  equations  for  the  ellipses  of  constant  likelihood  are  in  canonical  form)  as 


-vux+  ^-[(A  -  n)  ±  ((A  -  n)2  +  4i>2J]  w( 
f[}[(A  -  n)  ±  |(A  -  /i)2  +  4i>2J  J  +  *>' 


—  +  y  (A  —  n)  ±  |(A  —  ju.)2  +  4u2|  j  w, 

=  o  y/? 

Y“[(A  —  /x)  ±  |(A  —  (jl)'  +  4i/2|  J  +  V2 

where  the  choice  of  signs  has  the  same  effect  described  for  equations  (240)  and  (241). 
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